Cooling Dynamics of Photoexcited Carriers in Si Studied by Using Optical Pump and 
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We investigated the photoexcited carrier dynamics in Si by using optical pump and terahertz 
probe spectroscopy in an energy range between 2 meV and 25 meV. The formation dynamics of 
excitons from unbound e-h pairs was studied through the emergence of the ls-2p transition of 
excitons at 12 meV (3 THz). We revealed the thermalization mechanism of the photo-injected 
hot carriers (electrons and holes) in the low temperature lattice system by taking account of the 
interband and intraband scattering of carriers with acoustic and optical phonons. The overall cooling 
rate of electrons and holes was numerically calculated on the basis of a microscopic analysis of the 
phonon scattering processes, and the results well account for the experimentally observed carrier 
cooling dynamics. The long formation time of excitons in Si after the above-gap photoexcitation is 
reasonably accounted for by the thermalization process of photoexcited carriers. 
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I. INTRODUCTION 

Carrier dynamics in photoexcited semiconductors has 
been a central issue in semiconductor physics for 
decadesi. In particular, the carrier-phonon interactions 
have important effects in a wide range of optoelectronic 
devices, such as laser diodes and solar cells. Since the de- 
velopment of ultrafast lasers, numerous studies have been 
devoted to elucidating the photoexcited carrier dynamics 
by using nonhnear laser spectroscopic techniques such as 
optical pump and optical probe spectroscopy, four wave 
mixing, and transient gratings^. Time-resolved photo- 
luminescence has also been used to investigate hot car- 
rier relaxation and exciton formation^"— . With the devel- 
opment of mid-infrared ultrafast light sources, the wave- 
length of the optical probe has been extended to the lower 
energy region and it has become possible to access the 
intraband transitions, i.e., the free carrier responses, or 
the intra-exciton transitions with sub-ps temporal reso- 
lution. As such, the mid-infrared probe method has been 
used to study, e.g., electron- hole (e-h) droplet formation 
in CuCl^ and conversion dynamics of orthoexcitons to 
paraexctions in Cu20^. 

The advent of terahertz time-domain spectroscopy 
(THz-TDS) has enabled researchers to investigate the 
lower end of the terahertz frequency range — ^— . Un- 
like optical methods, THz-TDS enables the determina- 
tion of the complex optical conductivity from which one 
can quantitatively evaluate the carrier density and the 
scattering time of the photoexcited carriersi^"— , with- 
out resorting to the Kramers-Kronig relation. The tech- 
nique has been used to investigate the buildup of plasma 
screening in a high-density e-h system in GaAsi^ and 
formation of excitons in GaAs quantum wellai^^. 

Carrier dynamics in indirect gap semiconductors, 
in particular in Si, has also been extensively stud- 
ied. Near-infrared and visible pump-probe studies^^^— , 
and more recently time-resolved photoemission spec- 
troscopy^i^ have revealed the ultrafast dynamics of hot 
carrier relaxation at very high carrier densities (^J 10^^ 



cm^'^). However, the photoexcited dynamics at lower 
carrier density in Si is not fully understood yet, de- 
spite their inevitable importance in optoelectronic ap- 
plications. From the viewpoint of many-body quantum 
phenomena in e-h systems, an unambiguous understand- 
ing of the cooling dynamics of the photoexcited e-h sys- 
tem will be crucial for the realization of theoretically pre- 
dicted low-temperature quantum degenerate phases such 
as exciton Bose-Einstein condensation and the e-h BCS 
stated"—. A quantitative description on the slow forma- 
tion time of excitons on the order of several hundreds of 
ps, as observed in our previous study2^, is also needed 
for gaining an understanding of the fundamental optical 
properties of Si. 

The current study addressed the carrier thermalization 
dynamics of photoexcited e-h system in Si and the slow 
formation dynamics of excitons by using optical pump 
and THz probe measurements. We evaluated the time- 
dependent fraction of free carriers and excitons from the 
spectral shape analysis of the complex optical conductiv- 
ity in the terahertz frequency range and estimated the 
effective carrier temperature by using the Saha equation. 
This analysis was partly motivated by the scheme re- 
ported by Kaindl et al.^^. The cooling dynamics of the 
carrier system was analyzed numerically by calculating 
all the relevant phonon relaxation rates. The observed 
cooling rate was quantitatively reproduced without ad- 
justable parameters. Accordingly, the observed long for- 
mation time of excitons in Si after above-gap photoex- 
citation was well accounted for by the carrier thermal- 
ization dynamics dominated by electron (hole) - phonon 
scattering. 



II. EXPERIMENTS 
A. Methods 

We used the optical pump and THz probe technique 
to study the transient photoexcited free carrier responses 
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FIG. 1: (color online). Photo-induced change of dielectric function (left) and conductivity (right) at lattice temperature Tl (a) 
30 K and (b) 60 K for the excitation density of Ip—VJ fij/cm^ and for the indicated pump-probe delays At. The solid lines in 
each panel show the Drude model fitted to the low-energy region. 



and to investigate the formation dynamics of excitons 
through the observation of the ls-2p transition around 12 
meV (~3 THz)2£. The output from a multi-pass ampli- 
fier seeded by a 14-fs mode-locked Ti:sapphire oscillator 
was used as a light source. The duration of the amplified 
pulse was 30 fs with the repetition rate of 1 kHz, and 
the center wavelength was 800 nm (1.55 eV). The output 
from the amplifier was divided into three beams, respec- 
tively for THz generation, THz detection, and the optical 
pump. The THz probe pulse was generated from a 300 
fim thick GaP(llO) crystal by optical rectification of the 
laser pulse. The waveforms of the THz pulse transmitted 
after the sample with and without the optical pump were 
recorded by electro-optic sampling with a 300 /im GaP 
crystal2I. The complex transmittance change was then 
obtained by Fourier transformation, thereby yielding a 
photo-induced change in the complex dielectric function, 
Ae = Aei -I- iAai/eoUj, where Acti is the change in the 
real part of the optical conductivity. The frequency range 
of our experimental setup was between 0.5 THz (2 meV) 
to 6 THz (25 meV) . To make a homogeneous optical exci- 
tation in the lateral direction, the sample was irradiated 
with a pump pulse with a spot size of 3 mm in diameter, 
which was larger than the THz-probe spot size of d=4A 
(A ranged from 0.05 mm (6 THz) to 0.6 mm (0.5 THz)) 
in our experimental setup. 

As a sample, we used a high purity and high resistivity 
(10 kQ cm at room temperature) Si single crystal with an 
optically fiat (001) surface. To achieve a nearly homoge- 
neous optical excitation in the depth direction, the sam- 
ple was mechanically polished to a thickness of 22 fim, 
which was smaller than the penetration depth of 27 fim at 
1.5 eV. The sample surface was cleaned with hydrofluoric 
acid (HF) to remove the surface oxide layer. The same re- 
sults were also obtained without the HF treatment. The 



sample was freely mounted inside a cryostat. 

B. Formation dynamics of excitons 

We first studied the formation dynamics of excitons 
from the photoexcited unbound e-h pairs. The dots in 
Fig. [T]show the experimental results obtained at lattice 
temperatures of Ti=30 K (a) and 60 K (b). The exci- 
tation density was /p=17 /LtJ/cm^ and the delay time At 
between pump and probe pulses was varied. From the ex- 
citation pulse energy flux and the absorption coefficient, 
the excited e-h pair density N can be estimated to be (1.3 
± 0.4) X 10^^ cm""^. Please note that this value is be- 
low the Mott density at 30 K, Nmou = 7.4 x 10^^ cm^^^ 
calculated from the random phase approximation^^. The 
dielectric function and the conductivity spectra 10 ps af- 
ter the photoexcitation are well reproduced by the Drude 
model, according to which the dielectric function is ex- 
pressed as 



where eb = 11.7 is the background dielectric constant, 
7 the damping constant, Ng^ the free e-h pair density, 
eo the vacuum permittivity and to*=0.16too the optical 
mass of free carriers in Si-^, respectively. The Drude 
component decreases with increasing delay time, and a 
peak emerges around 10-12 meV, an energy that corre- 
sponds to the ls-2p transition of excitons in Si. To es- 
timate the free e-h pair density at each delay time, we 
can fit the data with the Drude model in the low energy 
region (2.0-7.8 meV), where the Drude component dom- 
inates the spectrum. The solid lines in Fig. [T] plot the 
fitting curves, without the contribution from the exciton 
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FIG. 2: (color online). The temporal change of the free e-h 
pair density Neh obtained from the fitting in Fig. [T] 



component. Figure [H shows the temporal change of N^h 
obtained by the Drude fitting. At Tl — 30 K, N^h con- 
tinues to decrease until about 400 ps. Contrary to this 
free carrier behavior, the exciton density, as indicated by 
the ls-2p absorption at 12 meV, gradually increases af- 
ter the photoexcitation and reaches a constant value after 
400 ps (Fig. [11(a)). Accordingly, the temporal decrease 
in Neh can be attributed to the formation of excitons. 
It should be noted here that the lifetime of the indirect 
excitons in Si is 1 /xs2!l, so that the dynamics observed 
in the current temporal region is free from the lifetime 
effect. Since it is considered that thermalization of the 
electron system occurs within 1 ps through the electron- 
electron collisions^, the observed long formation time 
of excitons can be presumably attributed to the thermal 
relaxation process of carriers in the lattice system. As 
shown in Fig. [21 a shorter decay time of about 100 ps 
was observed at = 60 K; i.e., the thermalization time 
of the carrier-lattice system is shorter than at Tl = 30 K. 
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FIG. 3: (color online). Photo-induced change in dielectric 
function (left) and conductivity (right) for the excitation den- 
sity of 17 fij /cia^ at various lattice temperatures. The solid 
lines in each panel show the Drude model fitted to the low- 
energy region. 
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C. Thermal ionization of excitons 

Since the exciton formation is complete at 400 ps at 
Tl > 30 K, we can see that excitons and phonons should 
reach thermal equilibrium after 400 ps. To confirm the 
thermal equilibrium condition, we set the delay time 
to Ai = 2 ns and investigated the stability of excitons 
against the thermal ionization by varying the lattice tem- 
perature (Fig. [3]). The excitation density was the same 
as for the results in Fig. [1] Ip=17 /iJ/cm^. With increas- 
ing Ti, the spectrum clearly changes from a Lorentz-like 
shape (excitonic ls-2p intraband transition) to a Drude- 
like one (free carrier), indicating exciton ionization due 
to thermal activation. To estimate the free e-h pair den- 
sity quantitatively, we performed a fitting of the Drude 
model following the procedure described in the previous 
section. The dots in Fig. [Hare plots of the resulting N^h 
as a function of Tl ■ Now, let us compare our experimen- 
tal results with the predictions of the Saha equationl^. 



FIG. 4: (color online). The lattice temperature dependence 
of the free carrier density Neh- The dots show the experimen- 
tally obtained N^h ■ The solid line shows the prediction of the 
Saha equation. 



The Saha equation describes the ratio of the density of 
excitons Nx and free carriers after statistical equilibra- 
tion of their chemical potential in the Boltzmann limit. 
For a 3D case, it is 

where /i=0.123mo and i?o=14.4 meV are the exciton re- 
duced mass"^'^ and the binding energy^. For a given total 
e-h pair density N — N^h + Nx , Eq. ^ yields the rela- 
tionship between carrier temperature Tc and free carrier 
density Neh- The solid line in Fig. [His the N^h obtained 



4 



100 




200 300 400 

Time (ps) 

FIG. 5: (color online). The temporal change in carrier tem- 
perature Tcit). The dots show the experimental data ob- 
tained from the results in Fig. [2] and the Saha equation. The 
solid lines show the calculated thermodynamics considering 
the carrier-phonon interaction. 




FIG. 6: (color online). The interband scattering processes in 
Si. g scattering means scattering to the opposite band, and 
/ scattering means a transition to one of the remaining four 
bands. 



Electrons 



from Eq. ^ with a fixed density 7V=l.lxlO^^ cm^^ (as 
estimated from the excitation pulse energy flux) and it 
shows good agreement with the experimentally obtained 
value of Nf.h. This agreement ensures the validity of our 
Drude model analysis and also the assumption that our 
experimental condition is at the Boltzmann limit (high 
lattice temperature and low excited e-h density). 



D. Cooling dynamics 

Considering that the free carriers and excitons can 
reach equilibrium within the experimental temporal res- 
olution of THz-TDS, ~ 2 ps, we can use Saha equation 
to deduce the transient carrier temperature Tc{t) from 
the transient free carrier density N^hit). The dots in Fig. 
[5] represent Tc{t), converted from Fig. [5] using the Saha 
equation. The plot for the higher lattice temperature 
Tl—60 K shows a fast thermal relaxation of carriers (^ 
100 ps), while the plot for the lower temperature shows a 
relatively slow one 400 ps). The reason for this differ- 
ence can be attributed to the phonon types involved in 
the thermodynamics, as will be discussed in the following 
section. 



III. ANALYSIS OF THE COOLING 
MECHANISMS 



To describe the cooling mechanism of the electron sys- 
tem from the microscopic point of view, we calculate the 
thermal relaxation dynamics considering the interactions 
between carriers (electrons, heavy-holes and light-holes) 
and phonons. 



For electrons, we consider the scattering process of in- 
traband acoustic phonons and interband phonons. The 
energy-loss rate per electron in a Maxwellian distribution 
of temperature Te is given by'^^ 



dt 



7r2/3 

X (fcBTe)3/2 ( 1 



t;4 

n p 



Tl 



(3) 



dt 



inter 



ttTi p 



(4) 



The subscript "ac" and "inter" denote the intraband 
acoustic deformation potential and interband scattering, 
respectively, mi and rrit are the longitudinal and trans- 
verse effective masses of the electrons^^, Sq is the intra- 
band acoustic deformation potential, and (DtK)^ is the 
interband deformation potential involving the phonon 
type z, the density of the crystal p=2.33 g/cve? , and the 
number of possible final equivalent valleys Zi = 1 {g scat- 
tering) or 4 (/ scattering) (see Fig. B{xq,xI) is 



huji 



(5) 
(6) 
(7) 



where uji is the phonon frequency for the phonon type 
i, and Ki is a Bcsscl function of the second kind with 
imaginary argument. Table U lists the physical parame- 
ters used in the calculation. 
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Because we consider the distribution function in the 
Boltzmann limit, is related to the average energy ac- 
cording to the equipartition theorem: 

{Ee) = y^BTe. (8) 

Using Eq. ([5]), we can calculate the cooling rate (dT^/dt). 
Figure [7] (a) shows the cooling rate per electron due to 
interband scattering, and Fig. [7] (b) shows a compar- 
ison of coolings due to intraband scattering and inter- 
band scatterings. As the phonon frequency increases, its 
contribution to the cooling rate becomes larger. In the 
present case, g scattering with LO phonons contributes 
the most to the cooling process. Figure [7] (b) indicates 
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FIG. 7: (color online). Calculated cooling rate per elec- 
tron, — dTe/dt, for a thermal (Maxwell-Boltzmann) distri- 
bution with electron temperature T^. The lattice temper- 
ature Tl is fixed to 30 K. Cooling rates due to different 
interband scatterings are shown in (a), where each line de- 
notes the following processes; g-LO: g scattering with longi- 
tudinal optical phonons, f-LA: / scattering with longitudinal 
acoustic phonons, f-LO: / scattering with longitudinal optical 
phonons, and Inter: the sum of above three contributions. A 
comparison of cooling rates between interband scattering and 
intraband acoustic phonon scattering is shown in (b), where 
each line denotes the following contribution; Inter: the sum 
of interband scatterings (corresponding to Inter in (a)). In- 
tra: intraband deformation scattering by acoustic phonons, 
and Total: the sum of all scatterings. 



TABLE I: 


Band parameters 


of electrons^ 




Value 


Units 


mi 


0.9163 


mo 


nit 


0.1905 


mo 


So 


9.0 


eV 




720 


IS. 


(DtK)g_LO 


11.0 


10**cV/cin 




550 


K 


(ntK)/_LA 


2.0 


10'*cV/cm 


(huj)f -TO 


685 


K 


(DtK)f-TO 


2.0 


10**cV/cm 



that the contribution from the intraband scattering be- 
comes larger than the interband scattering below 70 K, 
and the total cooling rate becomes smaller as de- 
creases. 

The present calculation does not include g scattering 
by TA and LA phonons or / scattering by TA phonons. 
Although these scatterings are forbidden in the selection 
rules^i^, there are several experimental indications^ 
that they are active in electron transport phenomena in 
Si. These claims are based on the fact that the selec- 
tion rules are evaluated for electron transitions between 
points exactly on the (100) axes, while, in general, elec- 
trons undergo transitions between points relatively far 
from these axes. In the present case, however, the dis- 
tribution function of electrons is not considered to be far 
from these axes because the density and temperature are 
relatively low, so we did not include these scatterings in 
our calculations. 



B. Holes 

Next, we consider phonon-scattering of holes. Because 
of the complexity of the valence band structure with its 
peculiarities of double degeneracy (heavy- and light-hole 
bands), nonparabolicity and anisotropy, the full descrip- 
tion of the transitions could be quite complicated. How- 
ever, for the present experiment, we can simplify the 
model as follows. Firstly, we consider only the intraband 
transitions for holes. Secondly, we include the acoustic- 
deformation (ac) and the nonpolar-optical-deformation 
(nop) potential as the scattering potentials. Thirdly, we 
replace the nonparabolic and anisotropic band structures 
with parabolic and isotropic ones with effective masses 
nii (i=light hole (Ih), heavy hole (hh)). (The validity of 
this simplification will be discussed later.) Accordingly, 
the energy-loss rates per Ih and hh in a Maxwellian dis- 
tribution of temperature Ti (i=lh, hh) are^ 

dE^\ 8{2)^/' Efm'/' 

dt ^2/3 ;,4^ 

x(fcsT.)^/^(l-^), (9) 
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dt 



nop 



2V'^ [DtKfm 



3/2 



TT/i p 



(10) 



for each scattering potential, where Ei and (DtK) are 
the acoustic- and the optical-deformation potential, re- 
spectively. Table |TT] lists the physical parameters used in 
the calculation. 



TABLE II: Band parameters of holes 





Value 


Units 


mih 


0.154 


mo 




0.523 


mo 


El 


5.0 


eV 


(?iLj)„op 


735 


K 


(DtK) 




10**cV/cm 



By applying Eq. ([8|) to the case of holes, we can also 
calculate the cooling rate {dTi/dt). Figure |5] shows the 
cooling rate per light hole (a) and heavy hole (b) in the 
form of a comparison of intraband and interband scatter- 
ing. The cooling rates behave similarly, as in the case of 
electrons. By comparing Fig. [5] (a) and (b), one can see 
that the heavy holes cool faster than light holes. This 
is because the cooling rate depends on the mass (Eq. 

Now let us discuss the validity of the simplified model. 
First, the nonparabolicity and the anisotropy of the va- 
lence bands should be taken into account in the high 
density or high temperature condition because the de- 
viation from the spherical and parabolic bands becomes 
pronounced when holes reach energies far from the bot- 
tom of the valence bands. Under our conditions, however, 
the density is not so high such that the chemical potential 
for holes is —9 meV with a carrier density of n = 1 x 10^^ 
cm"3 at 30 K. The thermal energy is fepT ~ 3 meV, and 
thus the energy of holes eh is distributed almost at the 
bottom of the valence bands compared to the spin orbit 
energy ego (eso=44 meV in Si). Therefore, the condi- 
tion Eh ^ Eso is fulfilled, and the valence bands can be 
approximated to be spherical and parabolic. Next, the 
treatment of interband and intraband transition in holes 
is accomplished by considering the overlap factor G{v), 
which is given by^ 



G{v) 



\ (l-t-3cos2 v) 



intraband transition 



I sin^ v : interband transition. 



(11) 



where v is the scattering angle. The interband transition 
that corresponds to the large-angle scattering can be ne- 
glected from the following reason. The large-angle scat- 
tering requires the participation of phonons having large 
wave number. Due to the linear dispersion relation, these 




FIG. 8: (color online). Calculated cooling rate per light hole 
(a) and heavy hole (b), —dTi/dt, for a thermal (Maxwell- 
Boltzmann) distribution with light and heavy hole tempera- 
tures Ti. The lattice temperature Tl is fixed to 30 K. Each 
line denotes the following scatterings; NOP: non-polar-optical 
phonon scattering, AC: acoustic phonon scattering. Total: the 
sum of these two scatterings. 



acoustic phonons with large wave number have high en- 
ergy. When Tl becomes low, the number of such high en- 
ergy acoustic phonons becomes small. Therefore, in our 
experimental condition (relatively small Ti,), the small- 
angle scattering is considered as the dominant process in 
the cooling mechanism. 



C. Cooling dynamics 

Having calculated each phonon-scattering processs 
contribution to the cooling rate, we can now calculate 
the thermodynamics. For this, we used an algorithm 
that was developed to describe the photoexcited carrier 
dynamics in GaAs^. The algorithm proceeds as follows: 
In the first step, we determine the initial carrier tempera- 
ture Tc'(O) from the excitation photon energy hv = 1.55 
eV, band gap Eg — 1.17 eV, and the involved phonon 
energy for the indirect optical transition in Si, ?iw=0.05 
eV: 



, ^ . N hv — E„ — huj 
ksTciO) = ;^ . 



(12) 
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The resultant initial carrier temperature 7c (0) is 2000 
K. In the second step, we calculate each temperature, 
i.e., for electron (Tg) and holes (T^ with i=hh,lh), after a 
finite time interval (Ai): 



-kBT,{t + At) ^ -kBTcit) 



dt 



dt 



^ + At, (13) 



inter 



At) = -ksTcit) 




dt 



At 



(14) 



nop y 



In our calculation, At is fixed to be 0.05 ps. In the third 
step, we unify the temperatures into one carrier temper- 
ature. The temperatures after At are different because 
of the difference of energy loss rates of electrons, light 
holes, and heavy holes. By taking into account that the 
scattering between carriers is much faster than the scat- 
tering between carriers and phonons, we must unify the 
temperatures at each step. The total energy density in 
each system is obtained by the average energy multiplied 
by its density {{E) x n). Considering that the average 
energy (Eq. ([5])) is proportional to the temperature, and 
that energy conserved in the carrier scatterings, the car- 
rier temperature is given by the following equation: 

The densities of the systems are related as follows: 



TVe = iVih + iVhh, 

iVih ^ /mihy/' 



(16) 
(17) 



Equation (I16|) is true from the charge neutrality condi- 
tion, and Eq. (IT71) is true by virtue of that in the Boltz- 
mann distribution, the density of states of ideal gas par- 
ticles having mass m in a 3D system is proportional to 
m^/"^. Thee fourth step returns to the second step with 
the new input parameter Tc{t + At) and the steps are 
repeated. 

The calculated results are shown by solid lines in Fig. 
[SJ and they are in very good agreement with the experi- 
mental data. Please note that no adjustable parameters 
were used. The rapid drop in temperature is due to inter- 
band scattering of electrons and optical phonon scatter- 
ing of holes. The subsequent slow decline in temperature 
is due to the intraband phonon-scattering of electrons 
and holes by acoustic phonons. The calculated cooling 
rates in Fig. [7] and Fig. [S] indicate that the cooling time 
of the carrier system becomes longer at lower T^. 



Finally, let us address the scattering processes ne- 
glected by our calculation. Scattering from impurities 
would give only a small contribution because we used 
high purity Si, from which only free excitons were ob- 
served in the terahertz spectra. The screening effect of 
electron-phonon interactions by the electron gas is known 
to be important in the high density regime^^, above 
Neh = 10^® cm~^. In our case, however, the density 
is much smaller, and thus we believe that the screen- 
ing effects are negligible. The simplified treatment using 
Maxwell-Boltzmann distribution function and also ne- 
glecting lattice heating may become inaccurate at very 
low Tl. Regarding this possibility, a more sophisticated 
calculation remains as future work. 



IV. SUMMARY 



We have investigated the cooling dynamics of photoex- 
cited carriers in Si by broadband terahertz time-domain 
spectroscopy. The formation dynamics of indirect exci- 
tons was revealed through the observation of the ls-2p 
transition of excitons around 3 THz. The relatively long 
formation time of excitons ranging from 100 ps to several 
hundreds of ps dependent on the lattice temperature is 
attributed to the thermal relaxation time of the carri- 
ers in the lattice system. To view the cooling dynamics 
quantitatively, we determined the carrier temperature as 
follows. First, we evaluated the fraction of free carri- 
ers and excitons from the spectral shape analysis of the 
complex optical conductivity. From the obtained ratio 
of the free carrier density to the total e-h pair density, 
we determined the transient carrier temperature by us- 
ing the Saha equation. To describe the observed cooling 
dynamics, we numerically calculated the carrier-phonon 
relaxation rate by taking into account all the relevant 
electron(hole)-phonon interactions in Si. We found that, 
at high temperatures, the carrier cooling process is dom- 
inated by the interband scattering of electrons with opti- 
cal and acoustic phonons and by the intraband scattering 
of holes with optical phonons. At low temperatures (be- 
low ~100 K), on the other hand, the contribution from 
the intraband acoustic phonons becomes dominant both 
for electrons and holes. The overall cooling rate of the e-h 
system was quantitatively calculated, and the calculated 
rate was in good agreement with the experimental result 
without adjustable parameters. This precise understand- 
ing of the electron-phonon interactions and the cooling 
mechanism of photoexcited e-h system in Si would be 
valuable for optoelectronics and photonics applications 
— , as well as for the realization of theoretically predicted 
low temperature quantum degenerate phases of e-h sys- 
temsSS. This work was partially supported by grants-in- 
aid (No. 22244036, No. 21104505), from MEXT, Japan 
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